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ABSTRACT 

We develop a model of gravitational lensing in a non-uniform plasma. When a grav- 
itating body is surrounded by a plasma, the lensing angle depends on the frequency 
of the electromagnetic wave, due to dispersion properties of plasma, in presence of a 
plasma inhomogeneity, and of a gravity. The second effect leads, even in a uniform 
plasma, to a difference of the gravitational photon deflection angle from the vacuum 
case, and to its dependence on the photon frequency. We take into account both effects, 
and derive the expression for the lensing angle in the case of a strongly nonuniform 
plasma in presence of the gravitation. Dependence of the lensing angle on the pho- 
ton frequency in a homogeneous plasma resembles the properties of a refractive prism 
spectrometer, which strongest action is for very long radiowaves. We discuss the ob- 
servational appearances of this effect for the gravitational lens with a Schwarzschild 
metric, surrounded by a uniform plasma. We obtain formulae for the lensing angle and 
the magnification factors in this case and discuss a possibility of observation of this 
effect by the planned VLBI space project Radioastron. We also consider models with 
a nonuniform plasma distribution. For different gravitational lens models we compare 
the corrections to the vacuum lensing due to the gravitational effect in plasma, and 
due to the plasma inhomogeneity. We have shown that the gravitational effect could 
be detected in the case of a hot gas in the gravitational field of a galaxy cluster. 

Key words: gravitation - gravitational lensing: strong - gravitational lensing: weak 
- gravitational lensing: micro - plasmas. 



1 INTRODUCTION 

The photon deflection angle in vacuum, in the Schwarzschild 
metric with a given mass M, is determined, for small deflec- 
tion angles a <C 1, by a formula 



4GM 
c 2 b 



2R S 
b ' 



(1) 



is a dispersive me dium with a permittivity tensor depend- 
ing on its density l|Landau fe Lifshit jfl960 ). In the disper- 
sive nonuniform medium a photon trajectory depends also 
on the photon frequency, and this effect has no relation to 
the gravity. The photon deflection in a non-homogeneous 
plasma, in presen ce of g ra vity, has been considered by 



where b is the impact parameter, and b ^> Rs, 
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is the Schwarzschild radius 
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1992). When the impact 



parameter is close to its critical value, corresponding to 
the capture of the photon by the black hole, the expression 
for the deflection a n gle is more complicated (iDarwinlfigHs " 



Misner et ail Il973l: Tvirbhadra fc Ellis! l2000l : iBozza et al 



200ll ; lBisnovatvi-Kogan fc Tsupkdl200a ) 



It is interesting to consider the gravitational lensing in 
a plasma, because in space the light rays mostly propagate 
through this medium. In an inhomogeneous plasma pho- 
tons move along the curved trajectories, because a plasma 
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in a linear approximation, when the two effects: the vacuum 
deflection due to the gravitation, and the deflection due to 
the non-homogeneity of the medium, have been considered 
separately. The first effect is achromatic, the second one de- 
pends on the photon frequency if the medium is dispersive, 
but equals to zero if the medium is homogeneous. 

A general theory of the geometrical optics in a curved 
space- t ime, i n arbitrary medium, is presented in the book of 
ISvngel (|l960l ). On the basis of his general approach we have 
developed the mode l of gravitational lensing in plasma . In 
our previous work ( Bis novatvi-Kogan fc Tsupkol |2009| ) we 
have shown, that in a nonlinear approach a new effect ap- 
pears. Even in the homogeneous plasma, the gravitational 
deflection angle really depends on the frequency of the pho- 
ton, and is different from the vacuum case. We have derived 
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the expression for the deflection angle of the photon in a 
weak gravitational field, in a weakly inhomogeneous plasma. 

In this work we use more general approach and derive 
the deflection angle for the photon moving in a weak grav- 
itational field, in the Schwarzschild metric, in the arbitrary 
inhomogeneous plasma. Such approach is more appropriate 
for the propagation in a cosmic plasma, because plasma den- 
sity changes significantly, from the density of the interstel- 
lar medium to the density of a black hole neighborhood. 
We consider here only the situation, when the whole deflec- 
tion angle, from the combined plasma and gravity effects, 
remains small. In the section 2 we derive a general expres- 
sion for the deflection angle in the inhomogeneous plasma, in 
the curved space-time. In the section 3 we discuss in details 
the important case of the weak Schwarzschild field with a 
spherically symmetric distribution of plasma, and derive the 
formula for the deflection in this situation. Our approach al- 
lows us consider two effects simultaneously: the difference of 
the gravitational deflection in plasma from the vacuum case 
and the non-relativistic effect (refraction) connected with 
the pl asma inhomogeneity. In the paper of iKulsrud fc Loebl 
( 1992) it was shown that in the homogeneous plasma the 
photon wave packet moves like a particle with a velocity 
equal to the group velocity of the wave packet, and with 
a mass equal to the plasma frequency. In the section 4 we 
show, that our result for a homogeneous plasma follows also 
from this approach. In the section 5 we discuss the obser- 
vational appearances for a Schwarzschild point-mass lens, 
surrounded by a uniform plasma (effect of gravitational ra- 
diospectrometer). In particular, we obtain formulae for the 
magnification factors in this case. We also estimate a pos- 
sibility of observation of this effect by the planned VLBI 
space project Radioastron. In the section 6 we consider the 
models with the nonuniform plasma distribution. For differ- 
ent gravitational lens models we compare the corrections to 
the vacuum lensing due to the gravity effect in plasma and 
due to the plasma inhomogeneity. 

After the publication of our previous work 
l|Bisnovatvi-Kogan fc Tsupkol l2009h two papers con- 
cerning the effect of amplification of gra vitational deflection 
in dispersiv e medium were published l|Dressel et al.l 120091 ; 
lLinetll2009l) . These papers have nothing to do with as- 
trophysics, authors calculate the vertical deflection of a 
light ray in a medium with strong frequency-dependent 
dispersion in a uniform gravitational field, and estimate 
such effect for the labor atory exper iments with the Earth 
gravity. In the paper of iLinetJ (|2009l ) the Synge method of 
calculation is used. 



2 DEFLECTION ANGLE IN THE 
INHOMOGENEOUS PLASMA IN 
PRESENCE OF GRAVITY 

Let us consider a static space-time with a metric 



o\2 



ds 2 — gik dx l dx k — g a p dx a dx^ + goo (cfe°) 
£,fc = 0,l,2,3, a,@= 1,2,3. 



(2) 



Here gik does not depend on time. Let us assume that the 
gravitational field is weak, so we can write 

gik = r)ik + hik, h ik -C 1, h ik — > under x a — > oo . (3) 



Here r/ik is the flat space metric (—1,1,1,1), and hik is a 
small perturbation. Note (|Landau fc Lifshitzll 19931 ). that 

ik ik iik ik iik 7 / i\ 

g =7) - h , rj =rj, k , h = h ik . (4) 

Let us consider, in this gravitational field, a static inhomo- 
geneous plasma with a refraction index n, which depends on 
the space location x a , and the photon frequency oj(x a ): 

' • Jj = *** N ^ =K B N{x*). (5) 



Here cj(x a ) is the frequency of the photon, which depends 
on the space coordinates x 1 , x 2 , x 3 due to the presence of 
the gravitational field (gravitational red shift). We denote 
w(oo) = uj, e is the charge of the electron, m is the electron 
mass, w e is the electron plasma frequency, N(x a ) is the elec- 
tron concentration in inhomogeneous plasma, and we do not 
assume that N(oo) = 0. 

The opt ics in a curv ed space-time, in a medium, was 
developed bv lSvngel (|l960l ). It was shown that, for the static 
case, the connection between the phase velocity u, and a 4- 
vector of the photon momentum p % , using of the refraction 
index of medium n, n — c/u, is written as 



1 + 



PiP 



(p V=m) 



(6) 



Here c is the light velocity in a vacuum. The refraction index 
n, defined for a plasma in ©, is a function of x a and uj(x a ). 
In the vacuum n = 1, and we can obtain from ([§} the usual 
relation for the square of the photon 4- vector: pip 1 = 0. In 
the medium, the square of the photon 4-vector is not equal 
to zero. For the medium in a flat space-time we have 



goo = -1, g aa = 1, p = -po, p a = ; 
-(p°) 2 + (p a ) 2 



1 + 



(P ) 2 



and obtain the usual relation between the space 
and time compon e nts o f the 4-vector of t h e photon 
llLandau fc Lifsliitzl 1 196(1 IZhelezniakovl Il977l ; iGinzburd 

liazeD, 



/ a\2 1, 0\2 

(P ) =n (p ) 



(8) 



For a static medium in a static gr avitational field, we have 
l|Svngelll960l . see also lMoller|[l972T ): 



Po 



(9) 



where h is the Planck constant. A zero component 
of the 4-momentum i s the energy divided by c (see 
lLandau fc Lifshitj[l993l ). so in a flat space-time we have: 



o 1 K 

Po = -p = rwj , 

c 



(10) 



where ui = w(oo). The components of the 4-vector p 1 , during 
an arbitrary motion in a non-homogeneous medium, in a flat 
space are written as 

/ Too nhco 



(p°,p1 



Pi = 



(po,p a ) = (- 



Tito nhuj 



(11) 
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where e a — p a /p and e a = p a /p (p = \/ p\ + p\ + P 2 ) are 
the unit 3- vectors in the direction of the 3- vector p a and p a 
correspondingly. In the flat space we have e a — e a - We see 
that for the photons moving in the plasma pip 1 = —m 2 ffC 2 , 
with m e ff = ^ y/1 — n 2 . Using (J6]) and ([9]), we see that this 
relation is valid also for any static gravitational field, with 
m e ff — ft "^| - Vl ~ n- 2 - We have then, using ((5J, that the 
effective photon mass m e f f and effective velo city v e ff in a 
plasma are written as l|Kulsrud fc Loebl (jl992T l) 

?uo(x a ) r huj e 

m -ff = -^vi - ™ 2 = — , 



2 2 

2 P C 

V eff - 



p 2 + m 2 ff c 2 



= 1 



uj 2 (x a ) 



2 2 2 

c — n c . 



(12) 



Thus the effective photon velocity equals to the group ve- 
locity of the photon in plasma 



d(nu)) 
duj 



(13) 



The energy of the photon in plasma E e ff, and relation be- 
tween the energy and the effective mo mentum p B ff are a lso 
the same, as for the massive particle ijMisner et al.lll973T ) 



E ef f = hui(x a ) = V— ffoo cp°, 



(14) 



Peff 



E, 



e f f 9 9 

2 - m e//C 2 



nhui(x a ) 



The trajectories of photons, in presence of the gravit ational 
field, may be obtained from the variational principle l)Svngd 
119601 ) 



Pi dx 1 I = 0, 



(15) 



with the restriction ([6]), which may be written in the form 



W{x\ Pi ) 



9 3 PiVj - {n 2 - 1) (po\/-g 00 ^ 



0. (16) 



Here we define the scalar function W(x , ',pi) of x % andpi. The 
variational principle (JTSj , with the restriction W(x l , P i) = 0, 
leads to the following system of differential equations l|Svngel 
Il960h : 



d^__m_ dp 1 _ dw 

d\ dpi ' d\ dx i 



(17) 



with the parameter A changing along the light trajectory. 
In the case of a plasma with the refraction index ([5|, the 
restriction ()16p can be reduced, with using of §5§ and @, to 
the form 



W(x\ Pi ) 



PiPj + 



o. 



(18) 



From (|17[) we obtain the system of equations for the space 
components x a , p a : 



dx a a p 
~dX= 9 p * 
or 
dx c 
1\ 

dp a 



dp a 

d\ 



■ 9,aPiPi 



If? 

2 c? 



(19) 



a/9 

= 9 Pp, 



,0, , 



1 h 2 Aire' 



dx -^P/^-^Po-^ — 



N, 



(20) 



It follows from (|17l) that in the static field the component 
po is constant along the trajectory. Let us consider a pho- 
ton moving along z-axis in a curved space-time, in an in- 
homogeneous plasma. Due to the curved space-time metric, 
and plasma inhomogeneity, the photon will move along the 
curved trajectory. We use the approximation, in which de- 
viations of the photon trajectory from the straight line are 
small. Therefore for the null approximation we use the 4- 
vector in a flat space 

i ( hw „ „ nhu)\ ( Hlj „ „ nhuj\ 

P \~' ' ' ~) ' Vl = \~' ' °'~) ■ ( ^ 

The unit 3-vector in the direction of the photon momentum 
is written in the null approximation as: e a = e a — (0, 0, 1). 
We integrate the equations (|20|l . calculating the right-hand 
side, by using the null approximation in the trajectory of 
the photon, with p % from (I21|l . We obtain than from the 
first equation in (|20|l 



dz nhui 
d\ c 



d\ = c 



dz 
nTiu) 



(22) 



The deflection angle is determined by a change of 3-vector 
e a , so let us express the second equation in (1 20 j) through the 
e Q . We obtain successively 



nhuj d(rihu e a 



dz 



1 I3j 

n 9,aPfiPl 



d(ne a ) dn de a 

dz a dz dz 



1 oo 2 1 fi w AT , ooN 
7:9,cPo- T:^K e N, a , (23) 



(24) 



de a 
dz 



dn 



" a dz 2 nh 2 uj 2 



h 

■ I gfaPpP-y + 9°"Po + -^K e N, a 



(25) 



In the right hand side of this equation we use the components 
from (|21|) . We are interested in the components of 3-vector 
e a , e3 = 1, which are orthogonal to the initial direction of 
the propagation, at a = 1,2. In the right hand side of ()25[1 
we use the null approximation with e a = 0. Using (Q, we 
obtain: 

de a 1 c 2 
~dz~ ~ 2n 2 h 2 u 2 X 



n 2 h 2 iu 2 h 2 u 2 

x I "33, a 5 h "00, a K~ 



h 2 



de a _ 1 
~dz~ _ 2 



f^33,a H K hoO,a 5 — ^K e N, a) 



(26) 



(27) 



And for the deflection angle a a = e a (+ao) — e a (— oo) we 
obtain: 



1 / / , hoo.a K e N >a , , 

a Q = ^ / I h 33 , a + r-Tj-jX ~ -T—2 I dz, 



1 — (cjI/oj 2 ) lo 2 — UJ 



a = 1,2. 



(28) 



For the deflection angle of the photon in the vacuum, iV = 0, 
we have 



= _ / ^ fl ' M ' a + ^oo,a) dz, a = l,2. 



(29) 
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For the axially symmetric problem, it is convenient to 
introduce the impact parameter b relative to the point 
mass, which remains constant in the null approximation 
for the photon moving along the axis z. The plasma has 
a spherically-symmetric distribution around the point mass, 
with the concentration N = N(r). In the axially symmet- 
ric situation the position of the photon is characterized 
by b a nd z, and th e absolute value of the radius-vector is 
r — yxf + x 2 + z 2 = s/b 2 + z 2 . We have the following ex- 
pression for the deflection angle in the plane perpendicular 
to direction of the unperturbed photon trajectory: 



1 

a b = - 



dr 



+ 



dhn 



K e dN(r) 



(coi/uj 2 ) dr uj 2 



dr 



dz. (30) 



Note that at < corresponds to bending of the light trajec- 
tory towards the gravitation center, and &b > corresponds 

to the opposite deflection. 

In our previous paper (jBisnovatvi-Kogan fc Tsupkol 

l2009h we have considered a weakly inhomogeneous plasma 
with N(x a ) = N + N^x 01 ), N = const, Ni < No- Here 
we assume that the deflection angle is small, but we do not 
assume that Ni is much smaller than No, so the representa- 
tion of N(x a ) as a sum is not required. 



3 LENSING IN THE SCHWARZSCHILD 
METRIC 

Let us calculate the deflection angle for a photon moving in 
the inhomogeneous plasma, in the Schwarzschild metric of 
the point mass M, with 



-c 2 (l-R s /r)dt 2 + -_ 



dr 2 



Rs/: 



■+r 2 (d8 2 +sin 2 8 dtp 2 ). (31) 



In the weak field approx imation this metric is written as 
jLandau fc Lifshit j|l993h 



ds 2 = dso + (c 2 dt 2 + dr 2 ), 
r 



(32) 



where dso is a flat part of the metric ds 2 = —c 2 dt 2 + dr 2 + 
r 2 (d8 2 + sin 2 8 dip 2 ). The components hn, are written in the 
Cartesian frame as ({Landau fc Lifshitzll 19931 ) 



Rs 



Rs 



Rs 2 a 
COS 8. 



(33) 



Here s a is a unit vector in the direction of the radius- 
vector r a — (xi,X2,%3), the components of which are equal 
to directional cosines, the angle 8 is the polar angle be- 
tween 3-vector r a =r a , and z-axis, and S3 = cos# = z/r = 
z/Vb 2 + z 2 . Using formula J30J) we obtain: 



Rs_ 
b 



Rsb 



(uj 2 /uj 2 ) r 3 



+ 



K e b dN(r) 



UJ — LUe r 



dr 



dz. 



(34) 



To demonstrate the physical meaning of different terms 
in (|34[) . we write this expression under condition 1 — n = 



uj 2 /lo 2 <C 1. Carrying out the expansion of terms with the 
plasma frequency, we obtain: 



ab = -- 



2R S 



1 Rsb 

2 uj 2 



1 Keb 

2 uj 2 



1 dN{r) 
r dr 



dz 



4dz- 



l K e b 
2 uj 4 



dN(r) 
dr 



dz. 



(35) 



The first term is a vacuum gravitational deflection. The 
second term is an additive correction to the gravitational 
deflection, due to the presence of the plasma. This term 
is present in the deflection angle both in the inhomoge- 
neous and in the homogeneous plasma, and depends on the 
photon frequency. The third term is a non-relativistic de- 
flection due to the plasma inhomogeneity (the refraction). 
This term depends on the frequency, but it is absent if the 
plasma is homogeneous. The forth term is a small additive 
correction to the third term. If we use the approximation 
1 — n — uj 2 /lo 2 <<iC 1, and neglect small second and the forth 
terms, we obtain a separate input of the two effects: the 
vacuum gravitational deflection, and the refraction deflec- 
tion in the inhomogeneous plasma. Calculation of the re- 
fraction deflection for a power-law concentration was give n 



in our previous work ([Bisnovatyi-Kogan fc Tsupkol 
see alsolBliokh fc Minakovl (119891)" "'" 



2009), 



1970), 



iMu hlcman c t al 
iThompson et all (|l994l '). iLightman et al.l i| 19791 ). Note that 
the refraction in the inhomogeneous plasma with N(r) — 
No(Ro I 'r) h , where iVo = const, Tio = const, h = const 7^ 0, 
leads to the refraction deflection angle a r , which is oppo- 
site to the gravitational deflection |Bliokh fc Minakovlll989l : 
iBisnovatvi-Kogan fc Tsupkol f2009). For u > oj e we have: 

_ 1 Ane 2 M /i?o\ h ^r(| + i) 



-No 



F(x) 



l dt. 



(36) 



For the arbitrary n one needs to use the expression (|34[) 
which is valid in a general case. The main approximation 
used here, is the smallness of the deflection angle, what can 
be satisfied even if the concentration iV changes significantly, 
or if the refraction index n is not close to the unity. The most 
interesting result following from our calculation is that even 
in the case of the homogeneous plasma the photon deflec- 
tion angle differs from the vacuum case, and depends on the 
plasma and photon frequency. Indeed, for uj e = const we 
obtain from (|34l) 



_Rs_ , 

'"' ~ b V 1-{u 2 /oj 2 ) 



(37) 



This formula is valid only for uj > lu b , because the waves 
with cj < uj e do not propagate in the plasma l|Ginzburd 
Il970h . Here at, < 0. This means that the light ray is bent 
in the direction of the gravitation center, as it occurs in the 
vacuum. Thus the presence of plasma increases the grav- 
itational deflection angle. This formula is valid under the 
condition of smallness of at, but this condition allows the 
second term in brackets to be much larger than the first one. 
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So, the gravitational deflection in plasma can be significantly 
larger than in the vacuum. This effect has a general relativis- 
tic nature, in combination with the dispersive properties of 
plasma. Such effect may happen only for radio frequency 
photons, because optical frequencies are much higher than 
the plasma frequency u) e , so the effect would be negligible. 

In the literature on gravitational lensing theory the de- 
flection angle is determined usually as a positive one fl]), and 
is defined as the difference between the initial and the final 
ray directions a = ej n — e ou t, where e is the unit tangent 
vector of a ray jSchneider et al.lll992h . Therefore, if we use 
this definition, we will have the expression with the opposite 
sign: 



_ Rs_ i , 

" b { + 1 - (wi/wa) I ' 



(38) 



which turns into the deflection angle for vacuum 2Rs/b, 
when uj — > oo. 



4 ANALOGY BETWEEN A PHOTON IN 
PLASMA, AND A MASSIVE PARTICLE 

Using the analogy between the motion of a photon in the 
homogeneous plasma, and a massive particle in the grav- 
itational field (|12p . (|14[1 we may easily find the deflection 
angle of the photon in plasma, produced by a point mass, in 
a weak field approximation. A test massive particle passing 
with the velocity v near a spherical body with a mass M, 
having the imp act parameter b S> Rs, deflects to the an gle 
a m defined as (Mi sner et al.|[l973l ; iLightman et al.lll979h : 



If we use in this formula v 



(39) 

cn = [l-( W ?y)] 1/2 



we obtain the same formula (|38|l for the photon deflection 
in the homogeneous plasma. Note that this analogy is valid 
for any static gravitational field, not only for a weak field. 



OBSERVATIONAL EFFECTS. THE 
PROJECT RADIOASTRON 



We see from (|38fl . that photons of smaller frequency, and 
larger wavelength, are deflected by a larger angle by the 
gravitating center. The effect of difference in the gravita- 
tional deflection angles is significant for longer wavelengths, 
when to is approaching oj e . That is possible only for the radio 
waves. Therefore, t he gravitational lens in plasma act s as a 
radiospectrometer l|Bisnovatvi-Kogan fc Tsupkoll2009h . 

The observational effect of the frequency dependence 
may be represented on the example with off-line lensing by 
the Schwarzschild point-mass lens, see fig.l, in the paper of 
i|Bisnovatvi-Kogan fc Tsupko| [2009) . This lens gives two im- 
ages of the source, on opposite sides of the lens. Angular po- 
sitions and forms of the images depend on the Schwarzschild 
radius of the lens, and on relative positions of the source, 
the lens, and the observer. The dependence of the deflection 
angle on the frequency in plasma leads to smearing of the 
images, and different parts of extended images have differ- 
ent spectra. In the case when the source, the lens and the 
observer are on the same straight line, the lensing image in 



vacuum is a thin circle l|Bisnovatvi-Kogan fc Tsupkoll2008h . 
In plasma this circle should have a finite thickness, with the 
radio spectra depending on the radius. The inner parts of the 
circle consist of a more energet ic photons than the its outer 
parts , see fig. 2, in the paper of (|Bisnovatvi-Kogan fc Tsupkol 
2009). In reality the gravitational lens is not a point mass, it 
has a complicated structure, and position of the images dif- 
fers from that of the point-mass lens. We should also notice 
that the source must be radio loud. 

The standard model of gravitational lensing is based on 
the Einstein deflection angle ([T}, which should be replaced, 
in presence of plasma by our formula (|38p . The angular half 
separation due to gravitati onal lensing, between the images 
of the source in vacuum (jSchneider et alJll992h is of the 
order of 



2R S 



D d 



D d D s 



(40) 



where Dd is the distance between the observer and the lens, 
D s is the distance between the observer and the source, Dds 
is the distance between the lens and the source. In the case 
of a perfect alignment of the source, the lens and of the ob- 
server, the image of the source is called Einstein ring, its 
radius has an angular size do- The observed angular separa- 
tion of quasar images is usually around 1 arcsec for lensing 
by a galaxy. Lensing in presence of a homogeneous plasma 
(|38[) leads to an angular half separation between images as 



apl 
7 o 




Rs 



D A 



DdD 3 



-"" w H 1+ i-(i/" 2 ) 



(41) 



which may be called, as plasma Einstein ring. For oj 2 /ui 2 <C 1 
we obtain 



npl 
y 



, 1 ul ^ 
1+ 4^' ""■ 



(42) 



The difference between angular separations of images in vac- 
uum and in plasma A#o, produced by the same lens config- 
uration, is equal to 



A0 O C - 



1 w£ 



2.0 ■ 10 7 



(43) 



where v is the photon frequency in Hz, w = 2-kv. The formula 
(|43|) gives the difference between the deviation angle of the 
radio wave with a frequency v, and the optical image, which 
may be described by the vacuum formula (|40|) . 

Let us estimate the possibility of observa- 
tion of this effect by the planned project Ra- 
dioAstron (see the RadioAstron web page at 
http:/ /www. asc.rssi.ru/radioastron/index. html ). The 
Radioastron is the VLBI space project led by the Astro 
Space Center of Lebedev Physical Institute in Moscow. 
The pay load is the Space Radio Telescope, based on 
the spacecraft Spektr-R. For the lowest frequency of the 
Radioastron, v = 327 • 10 6 Hz, the angular difference 
between the vacuum and the plasma images is about 10~ 5 
arcsec, when the plasma density on the photon trajectory, 
in the vicinity of the gravitational lens is of the order of 



6 G. S. Bisnovatyi-Kogan and 0. Yu. Tsupko 



5 ■ 10 4 cm 3 . this angular resolution is supposed to 



be reached in the project Radioastron. 

The magnification of the image increases with increas- 
ing of the deflection angle, therefore different images may 
have different spectra in the radio band, when the light prop- 
agates in regions with different plasma den sity. The magni- 
ficati on is determined by the deflection law l|Schneider et al.l 
1 19921 ). Let us demonstrate it on the example of the point- 
mass lensing. The magnification factor of the primary image 
fj,+ , located at the same side as the source relative to the lens, 
and of the secondary image fj,- , located at the opposite side, 
depend on the angular po sition of the sour c e. Co rresponding 
formulae can be found in Schneide r et all (1 19921 ) 



vV + 4 



Jy 2 + 4 
+ ^LE + 2 



^ = 4 



Vv 2 + 4 



Vv 2 + 4 



(44) 



(45) 



Here y = /3/9o, where /3 is the angular position of the source 
relative to the line passing through the observer and the 
lens. The total magnification of the source is equal to 



fltot = H+ + fi- 



y +2 



yvV + 4 



(46) 



Consideration of the total magnification factor is important 
for the microlensing events when the separated images are 
not resolved, and the only observable effect is changing of 
the flux from the source due to lensing. 

In the case of lensing in plasma we can rewrite these 
formulae, using y = fi/Q^ instead y: 



vi 



F- = 



v \Jv 2 + 4 

+ _v_y_ h2 



vV + 4 



vV+4 



- 2 



pi P i , P i y 2 + 2 
Kot = F+ + F- = 



(47) 

(48) 
(49) 



where y 



JL = £ 
e p Q l do 



1 / l 

2 I + 1 - (lj 2 /uj 2 ) 



-, -1/2 



V 



1 f 1 

2 I + 1 - (uj 2 /uj 2 ) 



-1/2 



(50) 



At large /3 the total amplification factor goes to unity, be- 
cause the influence of the lens on the light propagation be- 
comes negligible. For the Schwarzschild lens at a small angle 
f3, the amplification is inversely proportional to the angle /3, 
and is proportional to the angular radius of the Einstein 
ring, which increases with decreasing frequency, approach- 
ing infinity at the plasma frequency. As /3 goes to zero, the 
amplification increases, formally unrestrictedly for the point 
source. 

Let us consider the magnification in plasma, when 
1. Carrying out the expansion of the expression 
49l. and (1501) . we obtain the additional terms 



2 / 2 
UJ e /LO 



which arise in the case of the plasma, as compared to the 
case of the vacuum: 




Figure 1. The ratio fi^/fj,+ of the magnification of the primary 
image in presence of plasma to the same value in vacuum (lower 
curve) , and the ratio //^ //i— of the magnification of the secondary 
image in presence of plasma to the same value in vacuum (upper 
curve). Curves are plotted for w = y/2u} e . 



pi , 
= M+ + 



M- + 



y(y 2 + 4) 3 / 2 a; 2 ' 

i <4 

y(y 2 + 4) 3 /2 cj 2 ' 

2 w, 



(51) 



(52) 



(53) 



Vtat-Vtot + y(y2+ 4)3/2 u 2- 

We see that the presence of a homogeneous plasma increases 
the magnification. 

The ratio of the magnification of the primary image in 
presence of plasma n v l to the same value in vacuum /i+ is 
given in FigQ] for oj = s/2oj e , according to (|44)l . (|47)l . The 
upper curve in Fig.l is the ratio of the magnification of the 
secondary image in presence of plasma yP_ to the same value 
in vacuum /j,-, according to (145[l . (|48[) . As y approaches zero, 
the ratio of the magnifications f/£/fi+ and fi^/fi- becomes 
constant V+, depending on the frequency 



pi pi api 




(54) 



This value is equal »/§ at oj = V2uj e . At large y the ratio 



is tending to unity, the ratio fjF^/fi- is tending to 



the constant V- 



V- 



1 + 



{uI/uj 2 



(55) 



This value is equal 9/4 at cj = \/2uj e . Note, that both /x_ 
and fi_ are tending to zero at this limit. 

The deflection angle in presence of the plasma is larger, 
than in the vacuum, so the amplification for lower frequen- 
cies is larger, and in this situation the image spectrum differs 
from the original spectrum of the source, having more inten- 
sive the low- frequency part. The light in two lensing images 
is propagating through different media with different plasma 
density. Therefore, different images of the same source may 
have different spectra in the radio band. 
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In the ideal case, when two images of the same source 
are formed by rays propagating through the uniform plasma 
with different concentrations, the spectra of two images 
should be different. In the achromatic lensing the ratio of 
fluxes should be the same at all frequencies /x^ p< = 

^opt = ^rad^ go that ^pt^opt = jLl ™ < i^™d_ presence 

of plasma leads to larger amplification at lower frequencies, 
so if the ratio of fluxes in optics and in lower radio band 
of two lensing images is different /i° pt //i° pt 7^ / u+* d //ill a ' i , it 
may be related to plasma influence, so that the image with 
larger relative radio flux propagates through the plasma with 
larger density. In Fig. 2 two cases are shown for the value 
(M+/M+)/(M— /a 4 — )i as a function of Lj/u> e , for different fixed 
y. In the first one the plasma density is two times larger for 
rays forming the main image, and in the second case the ra- 
tio of densities is opposite. While the vacuum value /i+ 
does not depend on the frequency, one of these two behaviors 
is expected for the ratio of fluxes of two images as a function 
of the frequency. In the considered ideal case this plot should 
give the information about the ratio of plasma densities, and 
when the image with lower relative radio flux is formed by 
vacuum lensing, it would be possible to obtain the absolute 
value of the plasma density through which propagate rays 
from the image with higher relative radio flux. 

From Fig. 2 we see that the spectral dependences can 
have different forms for different y: while the dashed line 
is less than unity, the solid line can be both larger and less 
than unity. In Fig.2a in case when the light rays correspond- 
ing to primary image go through the densier plasma, we see 
that for smaller frequencies the primary image is more mag- 
nificated, relative to the vacuum magnification, due to the 
plasma presence than the secondary image, see solid line. 
In opposite case, when the light rays corresponding to the 
secondary image pass through greater density, we see that 
for the small frequencies the secondary image is more mag- 
nificated, therefore the dashed line is less than unity. For 
bigger y both curves are less than unity. In Fig. 2b, c we see 
that with increasing of y the solid line goes to the region, 
which is lower than unity. In Fig. 3 we plot the same de- 
pendence for the case, when the plasma density is ten times 
larger for rays forming the main image, and the second curve 
corresponds to the opposite ratio in densities. In the case, 
when the difference between densities is larger (comparing 
with situation in Fig. 2), effects connected with plasma are 
greater. Here the solid line remains above unity for greater 
values of y than for the situation described in Fig. 2. We 
should not e, that for y ^ 1 the second ary image is very de- 
magnified (Schneider ct al. 1992, 2006), so in a real situation 
we will have behaviour like on Fig. 2a and Fig3a, and ratio of 
densities can be uniquely defined. More detail can be found 
in Appendix 1. 

In this section we have discussed in details new ob- 
servational effects, connected with the gravitational deflec- 
tion in a homogeneous plasma. In reality the distribution of 
plasma around the gravitating objects is non-homogeneous, 
and there is a deflection connected with the medium inho- 
mogeneity. For the plasma density profile decreasing with 
the distance from center, the refraction deflection has a sign 
opposite to the gravitational one. Therefore the refraction 
and gravitational deflection in the nonuniform plasma par- 
tially cancel each other, and the value and the sign of the 
resulting angle depend on the particular configuration of the 



gravitational lens. In the following section we consider situ- 
ations of lensing with the non-uniform plasma distribution. 



6 MODELS WITH NONUNIFORM PLASMA 
DISTRIBUTION 

6.1 Singular isothermal sphere 

Let us consider a simple model of singular isotherma l 
sphere l|Binnev fc Tremaind 1 19871 ; IChandrasekharj Il939h 
which is often used in th e lens modeling of galax- 
ies and clusters of gala x ies dScrmeider et al.l 19921 . 120061 : 
iBartelmann fc Schneider! l200ll ; IWambsgansd I1998T I. The 
density distribution is written as 



p(r) 



(56) 



2nGr 2 ' 

Here a v is a one-dimensional velocity dispersion (for stars in 
galaxies or for galaxies in clusters of galaxies) . The projected 
surface mass density of a singular isothermal sphere to the 
lens plane, perpendicular to the light ray, is equal to 



v ; 2Gb 



(57) 



The deflection due to a axisymmetric mass distribution with 
the impact parameter b is equal to the Einstein angle for 
the mass M(b), where M(b) is the projected mass enclosed 
by the circle of the radius b. In another wor ds it is the 
mass inside the cylinder with the radius b (Clark Il972l : 
iBliokh fc Minakovlll989l ; ISchneider et~allll992l . l2006l v So we 
for this case we should rewrite formula (|35l) substituting 
M(b): 



Ot-b = 



AGM(b) 2GM(b) b 



c 2 b 



+ 



+ 



K e b f 1 dN(r) 
dr 



dz + 



K e b 



r dr 



(58) 







= Qi + «2 + Q?3 + Q4- 



Here we rewrite the formula in such a way that the gravi- 
tational deflection is positive as it is usually determined in 
the gravitational lensing theory (compare with (|35[) ). The 
term &\ is the vacuum gravitational deflection, the term &2 
is the correction to the gravitational deflection, due to the 
presence of the plasma, the term &3 is the refraction deflec- 
tion due to the inhomogeneity of the plasma, the term 04 is 
a correction to the third term. We are interested mainly in 
the effects, described by the terms &i, ch and 0:3. 

For the singular isothermal sphere the projected mass 
M(b) is: 



b 



M(b) = / E(fo') 2-Kb'db' = j$ Tib. 



(59) 



So the vacuum gravitational deflection due to the singular 
isothermal sphere is constant: 

q 1 =4tt4- (60) 
For the concentration of the plasma we have: 
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a) 



^Ai + )%. p V) 



b) 



1.5 2 



co/co, 



cO 



C) 




(^+ P V + )%_ P V_) 



y = 0.8 



5 E 7 8 9 



3 4 5 6 7 8 9 

co/co e o 



a) 




(^ + % + )/( M . p 7n_) 



y = 0. 



3.5 4 5 6 7 3 9 10 11 12 13 14 15 20 25 

co/o) c o 



b) 



y = 3 



35 4 5 B 7 S 9 10 11 12 13 14 15 20 25 



C0/C0 e 



C) 



y=12 



3.5 4 5 6 7 



9 10 11 12 13 14 15 20 25 

C0/C0e0 



Figure 2. The value ), as a function of uj/u> e , 

for different fixed y. w e + is the plasma frequency corresponding to 
the plasma concentration JV_|_ , which rays forming primary image 
go through. iu e - is the plasma frequency corresponding to the 
plasma concentration 7V_ , which rays forming primary image go 
through. The solid curve is plotted for uj e + = cj e oV2, = t^eO 
(o; e o is a dimensional constant), it corresponds to N+/N— = 2. 
The dashed curve is plotted for o; e + = uj e o, u> e - = w e oV2, it 
corresponds to N+/N— = 1/2. 



Figure 3. The value (/i J ^/^t+)/(/i^ i //i_ ), as a function of uj/u> e , 
for different fixed y. The solid curve is plotted for u) e + = ui e o \/To, 
Lj e — = ui e o (uJeO is a dimensional constant), it corresponds to 
N+/N— = 10. The dashed curve is plotted for w e + = u> e o, Ue— = 
oJeQvW, it corresponds to N+/N— = 1/10. 



N(r) 



p(r) 



(61) 



where m p is the proton mass, and n is a non-dimensional 
coefficient responsible for the dark matter contribution, and 
is approximately equal k ~ 6. We assume here that p(r), 
which is given by the formula ([56} , is the density of all kind 
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of matter, not only plasma particles. Thus, the plasma fre- 
quency is equal to 



= K e N(r) = _i2-p( r ). 



(62) 



The correction to the gravitational deflection due to the 
presence of the plasma is: 

oc 

. _ b 2 a 2 v K e Mr), 



(63) 



L0 2 C 2 K,m p J r 3 




dz 



p G J (b 2 + z 2 ) 5 / 2 ' 



Integration of such express ions can be performed using 
iGradshtevn fc Rvzhikl (| 1965h . and the properties of the T- 
function: 



dz 



i ^r(f + l) 



where T(x) = / t x ~ x e dt. 
o 

Using l|64|l. we obtain the deflection angles &2 and ay. 



(64) 



2 at Ke 



Q-2 — 



a 3 = 



3 C 2 KTUp GbJ 2 b 2 ' 



Keb at 



u) 2 K,m p tyG J (z 2 + b 2 ) 



dz 



1 K e 



4 KTTlp Gu) 2 b 2 



For the ratio of the angles &2 and 0:3 we have: 

6 2 _ 8 al 

63 3 c 2 ' 



(65) 
(66) 

(67) 



6.2 Non-singular isothermal gas sphere 

Let us consider a gravitational lens model of an isothermal 
sphere, in w hich the singularity at the origin is replaced by 
a finite core (|Hinshaw fe Krausdll987l ; IWufll996h : 



p(r) = 



2nG(r 2 +r 2 ) 



' p °-2^!> (68) 



where r c is the core radius. The corresponding projected 
surface mass density for this model is 



E(6) = 



2GVb 2 + r 2 ' 



(69) 



The total projected mass within b and the vacuum gravita- 
tional deflection angle are: 



M(b) = *§(V- 



&l = 47T 



b 2 + r'i - r< 



al Vb 2 + r 2 - r c 



(70) 



(71) 



Analogically to the previous subsection, we reduce the angle 
&2 to the form: 



a.2 = 



2GM(b)b K a a 2 v 



dz 



Km v 2ttG J r 3 (r 2 + r 



(72) 



2GM(b)b K e al 



dz 



c 2 co 2 nm p 2-kG J {z 2 + b 2 ) 3 / 2 (z 2 + b 2 + r 2 ) ' 


By substitution z 2 j(z 2 + b 2 ) = x 2 the integral 

dz 



h 



I 





(z 2 + b 2 ) 3 / 2 {z 2 + b 2 + r 2 ) 



is reduced to 
1 

h - 



1 

6V 2 



dx. 



(73) 



rl{b 2 +r 2 (l - x 2 )) 
Integration can be performed with using of lDwightl (|l96lh : 







b 2 r 2 



arctanh 



Vb 2 + r'i 



For «2 we obtain 
at K 



012 



c 2 uj 2 G Km 



1 

feV 2 



^bU. 

Kirip V 



b 2 +'. 



(74) 



(75) 



arctanh 



Vb^T' 



Angle &3 can be calculated using formula 

00 

K e b 1 a 2 . f dz 



Q3 



uj 2 K,m p 2nG J (z 2 + b 2 + r 2 ) 2 



1 al Ke b 



(76) 



4cj 2 G nm p (b 2 + r 2 ) 3 / 2 ' 

Let us express angles &2 and Q3 in terms of the central 
density po, using (|68p . We obtain for r c 3> b: 



2ix 2 Gpl Ke , 
a2 = — 5—5 br c , 03 



vrpo K e b_ 
2lj 2 Km v r c 



(77) 



and for 6>r c , similar to the case of the singular isothermal 
sphere, we have 



8n 2 Gp 2 Ke r 4 c „ 

OL2 = - - T7T, a?3 = 



TTpO Ke rl 



3c 2 



.i - nm p b 2 



(78) 



Introducing the mass of the uniform core Mc = ±fp r 3 c , and 
its gravitational radius R sc = 2G * /c , we obtain the ratio of 
these angles as 



3R 3C 
~2rT 



{r c > b), 



2R S 



(r- c < b). 



(79) 



In the realistic cases ||jp-| <C 1 in (|67[) . (|79[) . because 
spheres have a v <C c, and R sc <C r c . Besides, these re- 
lati ons are needed for the stability of i sothermal spheres, 
see iBisnovatvi-Kogan fc Zeldovichl l|l969l ). Therefore in this 
configuration the nonuniform plasma deflection effects are 
much stronger, than the gravitational plasma effects, and 
have an opposite direction. 
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6.3 Plasma sphere around a black hole 



6.4 Plasma in a galaxy clusters 



Let us consider a black hole of mass Mo, surrounded by 
the electron-proton plasma. We will consider a case, when 
we can neglect the self-gravitation of the plasma particles, 
compared to the gravity of a central black hole. Let us find, 
in the newtonian approximation, a density distribution of 
the plasma in the gravitational field a central point mass 
Mo. The equation of hydrostatic equilibrium for spherically 
symmetric mass distribution of a isothermal gas with the 
eq uation of state P = P$tT i n the field of the cent ral mass 
is (jBinnev fc Tremaindll987l ; IChandrasekhaill 19391 ) 



KTdp _ GMp 
p dr r 2 



(80) 



Here p, P and T are the density, the pressure and the tem- 
perature of the plasma, 5R = fc_g/m p is the gas constant, fes 
is the Boltzmann's constant and m p is the proton mass. To 
obtain a deflection angles we need to find the plasma den- 
sity distribution from the equation (|80[) and calculate the 
integrals with this density in the formula (|58[) . But it is in- 
teresting that we can find the ratio of angles Icb/etal without 
solving the equation (|80[) and calculation of the integrals in 
&2 and <3f3. Let us rewrite the equation (|80[) to the form: 



K e b^ T ldN(r) 
ui 2 r dr 



K e b GM , 



(81) 



where N(r) = p(r)/m p is the plasma concentration. If we 
compare it with expressions for &2 and 03 in (|58[1 and take 



into account that ui e = (4-7re /m)N(r) = K e N(r), we obtg 



that 

«2 
&3 



25RT 



(82) 



We see that this ratio does not depend on the value of the 
central mass Mo. For the boundary condition 



p(ro) = Po 



(83) 



we obtain from the equation (|80l) thc density distribution 



omq ( 1 1 ' 

p(r) = p e !KT \ r r o, 



Poe 



»(*"*) 



,B=™ (84) 



5RT 



For angles cb and &3 we obtain: 
2GM 6 K e 



Q-2 — 



Q3 



■ po e T o Int2, 



K e b 



B poe r o Int2, 



where Int2 — 



(b 2 + Z 2 )3/ 2 



b +* dz. 



(85) 
(86) 

(87) 



Calculation of Int2 can be found in Appendix 2. 

While the temperature of the non-self-gravitating 
sphere may have arbitrary values, the plasma effects may be 
comparable, and even less than the GR plasma effects. It is 
due to the fact, that with increasing temperature the plasma 
density can become arbitrary uniform, with corresponding 
decreasing of non-uniform plasma effect for refraction. We 
have used newtonian non-relativistic description of the gas 
sphere, but from the arguments listed above it is clear, that 
this conclusion remains valid also in the correct relativistic 
consideration. 



In a galaxy cluster the electron distribution may be more 
homogeneous due to large temperature of electrons. An ap- 
propriate approach for this case is to consider a singular 
isothermal sphere as a model for the distribution of the grav- 
itating matter, neglecting the mass of plasma, and to find a 
plasma density distribution from the solution of the equation 
of the hydrostatic equilibrium of plasma in the gravitational 
field of a singular isothermal sphere. In this approximation 
the density distribution of the gravitating matter has a form 

Mr) - ^ (88) 
with the vacuum gravitational deflection angle 

a!=47r4- (89) 



For the mass inside a sphere with a radius r we have 



(90) 



and the equation of a hydrostatic equilibrium of the plasma 
is: 



(91) 



UT dp _ 2a 2 v 
p dr r 

For the boundary condition p(ro) = po we obtain the plasma 
density from the equation (|91[) as 

'« = S = W- ^ 

The angles &2 and 63 are calculated using formula (|64|) . so 
we have 



Q2 



Q3 



2a 2 v n^ K e T(f + 1) 



c 2 uj 2 (s+l)m p r(* 



■po 



(?)" 



K s b 1 



■ Po r s 



dz 



(z 2 +b 2 )"/ 



2+1 



(93) 



(94) 



K e y/iv r(§ + 1; 



■po 



LJ 2 m p T(f) 
The ratio of these angles is equal to 



r(f + i)F(f) 
r(f + |)r(f + |; 



(95) 

2 1 2 1 

Under the condition s< 1 what corresponds to 2a 2 <C 5RT, 
the expressions are simplified to 

,; - = 2-k-^— po [ ^ ) , (96) 



2 



p 

Ke 



«3 = 



#tTuj 2 m 1 
25RT 



(97) 



(98) 



If relativistic plasma is present in a galaxy cluster, for exam- 
ple in jets from AGNs, the plasma GR effects may be larger 
than the effects of the nonuniform plasma. If the distribu- 
tion of plasma is not spherically symmetric, there may be 
distribution of plasma with the density gradient opposite to 
the direction of the gravitational force, for example, in the 
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presence of rotation. In this situation the angles &2 and 0:3 
may be of the same sign. 

Let us estimate Thomson optical depth for Thomson 
scattering of photons in the case of gravitational lensing in 
plasma, and depth of the bremsstrahlung absorption. For 
Thomson optical depth r we have: 

r « Na e L, (99) 

where N is the electron concentration, a e — 6.65 • 10 -25 
cm 2 is Thomson scattering cross section, L is characteristic 
distance for scattering. So the optical depth will be r < 1 

for N < 0.5 ■ 10 4 cm" 3 . 

For optical d epth for the abs orption of radiation in case 
cj 3> oj e we have (|Ginzburelll970h : 

j.3/2" 



17.7 + In- 



T 3/2 V 2 

So the optical depth will be r < 1 for 



N ~ 0.2 ■ 10' 



Go 4 k) 



3/4 



10 9 Hz 



10 2 pc 



1/2 



(100) 



(101) 



7 CONCLUSIONS 

In an observation of two images of the lensing point source in 
presence of plasma, the following features may be observed. 

1. Spectra of two images may be different in the long 
wave side due to different plasma properties along the tra- 
jectories of light rays forming the images. 

2. The extended image may have different spectra in 
different parts of the image, with a maximum of the spec- 
trum shifting to the long wave side in the regions with a 
larger deflection angle. 

3. The presence of plasma may influence on the timing 
effects in binary relativistic systems, similar to the double 
pulsar system J0737-3039A - J0737-3039B, and may induce 
the spectral dependence of the properties of the fluctuations 
of the microwave background radiation. 

4. We have carried out the calculations for models 
with the nonuniform plasma distribution: singular and non- 
singular isothermal sphere; for hot gas inside the gravita- 
tional field of a black hole, and of a cluster of galaxies. 

5. For different gravitational lens models we compare 
the corrections to the vacuum lensing due to the gravity 
effect in plasma, and due to the plasma inhomogeneity. We 
have shown that the gravitational effect could be detected 
in the case of a hot gas in the gravitational field of a galaxy 
cluster. 

6. We made estimations of the optical depth due to 
Thomson scattering and free-free absorption in the process 
of during the gravitational lensing in plasma. 



APPENDIX 1. INVESTIGATION OF 
SPECTRAL BEHAVIOUR IN FIG.2,3. 

Value (/i+ /a 4 -) under oj 3> 0J e has the following 
form: 



(m?//*+)/0£/a*-) = i + 



pi, 



2uj 2 (y 2 + 4) 



(102) 



(y 2 + 2 - y vV + 4) 



-to 2 e _ (y 2 + 2 + yv V+l 



We see that: 

a) if the rays forming the primary image go through 
greater density than the rays forming the secondary image, 
that is iu e + > (jj e -, then value (fjF / fj,+) / (fjF_ / fj,- ) can be 
both larger and less than unity (under big uf) depending on 
the angular position of the source y. In Figures 2,3 we see 
both situations, the solid lines) 

b) if the rays forming the secondary image go through 
greater density than the rays forming the primary image, 



that is cj e + < We-, then the value (//j. /fi+)/(jj, p _ //i-) can 
be only less than unity (see Figures 2,3, the dashed lines). 

Critical value of y, separating two cases of the location 
of the solid line (below unity and above unity) equals to 



V-8 + 6^/2 ^ nl<0 
ycrit = = 0.348 

for w 2 + = 2oj 2 _ (the solid lines on Fig. 2), and 



ycrit — 



y/-200 + llOyTo" 
10 



1.216 



(103) 



(104) 



for u> 2 + = Wlo 2 _ (the solid lines on Fig. 3). 



APPENDIX 2. CALCULATION OF THE 
INTEGRAL INT2. 



By substituting 1 / z 2 + b 2 = x the integral Int2 is reduced 
to 



i/b 



Int2 



x Bx , 

: e ax. 



V¥^ 2 



(105) 



Integration can b e performed with using of 

iGradshtevn fc Rvzhikl (|l965l ): 

Int2 = 1 (l + J [h(B/b) + Li{B/b)]j , (106) 

where 7i is the Bessel function of the first kind and L\ is 
the Struve function. 

Under condition B/& <C 1 which corresponds to 
GMo/b < KT, the integral Int2 is reduced to 
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